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An explicit model of F(R) gravity with realizing a crossing of the phantom divide 
is reconstructed. In particular, it is shown that the Big Rip singularity may appear in 
the reconstructed model of F(R) gravity. Such a Big Rip singularity could be avoided 
by adding R 2 term or non-singular viable F(R) theorjfflto the model because phantom 
behavior becomes transient. 
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1. Introduction 

It is obser vationally sup ported that the current expansion of the universe is 

1010 a c (i 7 o Tnl 

accelerating! | | | ,J | U | ' | 0|a | A number of scenarios to account for the current accel- 
erated expansion of the universe have been proposed (for reviews, see Refs. [TU1 [TT1 
[ia[T3[Tl[l5l[Ti[ni[T8l[T9l). 

Approaches to explain the current accelerated expansion of the universe fall 
into two broad categories. One is the introduction of some unknown matter, which 
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is called "dark energy" in the framework of general relativity. The other is the 
modification of the gravitational theory, e.g., "F(R) gravity", where F(R) is an 
arbitrary function of the scalar curvature R (for reviews, see Refs. [14] [15j [T6l [T7J 

□sua. 

Recent various observational 

dat£ J20|21|22|23| 

imply that the effective equation 
of state (EoS), which is the ratio of the effective pressure of the universe to the 
effective energy density of it, may evolve from larger than —1 (non-phantom phase) 
to less than —1 (phantom one), namely, cross —1 (the phantom divide). 

Various investigations to realize the crossing of the phantom divide have been 
executed in the framework of general relativity: Scalar-tensor theories with the non- 
minimal gravitational coupling between a scalar field and the scalar curvature or 
that between a scalar field and the Gauss-Bonnet term, one scalar field model with 
non-linear kinetic terms or a non-linear higher-derivative one, phantom coupled to 
dark matter with an appropriate coupling, the thermodynamical inhomogeneous 
dark energy model, multiple kinetic k-essence, multi-field models (two scalar fields 
model, "quintom" consisting of phantom and canonical scalar fields), and the de- 
scription of those models through the Parameterized Post-Friedmann approach, or 
a classical Dirac field or string-inspired models, non-local gravity, a model in loop 
quantum cosmology and a general consideration of the crossing of the phantom 
divide (for a detailed review, see Ref. [T3|) . However, explicit models of modified 
gravity realizing the crossing of the phantom divide have h ardly bee n examined, 
although there were suggestive and interesting related works CHIMES] 

In the present paper, we review our results in Ref.[5S]and reconstruct an explicit 
model of F(R) gravity in which a crossing of the phantom divide can be realized 
by using the reconstruction method proposed in Refs. [37] [25] (for more detailed 
references, see references in Refs. [55] [Ml 1301 (S3 l3"2"]) . It is demonstrated that the 
Big Rip singularity may appear in the reconstructed model of F(R) gravity. 



2. Reconstruction of a F(R) gravity model with realizing a 
crossing of the phantom divide 

2.1. Reconstruction method 

To begin with, we briefly review the reconstruction method of modified gravityPZEFJ 
The action of F(R) gravity with general matter is given by 

■F(R) 



S 



d Xy/—g 



2k 2 



(1) 



where g is the determinant of the metric tensor g^ v and £ m atter is the matter 
Lagrangian. 

By using proper functions P(<f>) and Q(4>) of a scalar field 0, the action in Eq. ([T]) 
can be rewritten to the following form: 

1 



5' 



2h- 



matter 



(2) 
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The scalar field <j) may be regarded as an auxiliary scalar field because <f> has no 
kinetic term. From the action in Eq. ([1]), the equation of motion of 4> is derived as 

dP(4>) dQ{4>) 

which may be solved with respect to <j> as (f> = <p{R). Substituting <j> = 4>{R) into 
the action ([2]), we find that the expression of F(R) in the action of F(R) gravity in 
Eq. (TTJ) is given by 

F(R) = P(<f>(R))R + Q^(R)). (4) 
From the action in Eq. ([2]), the gravitational field equation is given by 

\ 9la , [P(4>)R + Q(<P)] - R^P{4>) - g^ap(<p) + v„v„p(0) + k 3 ^*^ = o , (5) 

where V M is the covariant derivative operator associated with g^ v , □ = <? Miy V M V„ 
is the covariant d'Alembertian for a scalar field, and j^™ attcI ) j s contribution to 
the matter energy-momentum tensor. 

We assume the flat Friedmann-Robertson- Walker (FRW) space-time with the 
metric, 

ds 2 = -dt 2 + a 2 (t)dx 2 , (6) 

where a(t) is the scale factor. 

In the flat FRW background ((6]), the (/x, v) = (0,0) component and the trace 
part of the = component of Eq. ([5j), where i and j run from 1 to 3, 

become 

- 6H 2 P(q>(tj) - Q{cf>{t)) - 6H dPi j t (t)) + Up = , (7) 

and 

o d2p (m , AIT dP(<Kt)) 



dt 



2 



4H dt + + 6H2 ) P(0(<)) + O(0(<)) + 2K2P = ' (8) 



respectively, where H = a/a is the Hubble parameter and a dot denotes a time 
derivative, ' = d/dt. Here, p and p are the sum of the energy density and pressure 
of matters with a constant EoS parameter Wi, respectively, where i denotes some 
component of the matters. 

By eliminating Q((j>) from Eqs. ([7|) and ((8]), we obtain 

«_ ff « +2J j m())+ ^ +p) = „. (9) 

We note that the scalar field <f> ma y be taken as cj) = t because <j> can be redefined 
properly. 

We consider that a(t) is expressed as 

a(t)=aexp(g(t)) , (10) 
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where a is a constant and g(t) is a proper function. In this case, Eq. © is reduced 
to 

cPP{4>) dg(</>) dPtf) <Pg(<l>) 
d<P 2 d<j> d<t> d<t> 2 [<P) 

+ n 2 (1 + «"<) p t a-^ 1+w ^ exp [-3 (1 + Wi ) ~g{4>)\ = , (11) 

i 

where pi is a constant and we have used H = dg(cf>)/ (d<f>). Moreover, it follows from 
Eq. that Q(<f>) is given by 

2 



Q(4>)=-6 



dm 



P(4>) - 6 



dg(<f>) dP{4>) 



+ 2k 2 Pia~ 3(1+Wi) exp [-3 (1 + Wi ) g{cj>)] . (12) 

i 

Hence, if the solution of Eq. (fTTj) with respect to P{4>) is obtained, we can find Q{(j))- 



2.2. Explicit F(R) gravity model with realizing a crossing of the 
phantom divide 

Next, we reconstruct an explicit model of F(R) gravity in which a crossing of the 
phantom divide can be realized by using the reconstruction method explained in 
the preceding subsection. 

A solution of Eq. (ITU) without matter is given by 



P(0) = e ^ )/2 p(<£) 
g{<t>) = -10 In 



to 



7+1' 



P((f>)=P+(t> l3+ +P- 



l±y/l + 1007(7+1) 



(13) 

(14) 

(15) 
(16) 



where 7 and C are positive constants, to is the present time, and p± are arbitrary 
constants. 

From Eq. (Q~4|), we see that g(4>) diverges at finite when 

= ^ = t o C- 1 /^+ 1 ) , (17) 

which implies that there could be the Big Rip singularity at t = t s . 

We consider only the period < t < t s because g(4>) should be real number. 
From Eq. (|14[) . we obtain the following Hubble rate H{t): 

"7(^)"' 1 + (7 + l)gfe)" 



H(t) = 



dm 



7+1 



(18) 



where it is taken d> = t. 
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In the flat FRW background ©, even for F(R) gravity described by the action 
in Eq. ([1]), the effective energy-density and pressure of the universe are given by 
Pes = 3H 2 / k 2 and p c $ = — (^2H + 3H 2 ^j / k 2 , respectively. The effective EoS w, 

Pes/Pes is defined a 



"off 



W c ft 



2H 
3H 2 



(19) 



For the case of H(t) in Eq. (|T8l) . from Eq. (fT9| we find that w e s is expressed as 

w eS = -l + U(t), (20) 

where 

/ \ 2 7 + l / \ 

_ 7 + 4 7 ( 7 + l)(±J +( 7 + i)^j 



,2(2 7 +l) 



U(t) 



2H 
3H 2 



15 



7 +(7 + !)(£) 



\ 27+1 



(21) 



Furthermore, the scalar curvature is given by R = 6 (h + 2H 2 ^j . For the case of 
Eq. (fT8|) . R is described as 



GO 



R 



/ \ 27+1 / \ 

7 (20 7 - 1) + 44 7 (7 + 1) {J-J + (7 + 1) (20 7 + 21) {± J 



2(27+1) 



t 2 



(A) 



27+1 



(22) 



In deriving Eqs. ^2TJ) and ([22|l. we have used Eq. (fT7| . 
When t — > 0, i.e., t -C t S j ^(*) behaves as 

»( ( >~i?. 

In this limit, it follows from Eq. (TT91 that the effective EoS parameter is given by 

1 



(23) 



Weft 



-1 



157 



(24) 



This behavior is identical with that in the Einstein gravity with matter whose EoS 
is greater than — 1. 

On the other hand, when t — > t», we find 



H(t) 



10 

ts-t 



(25) 



In this case, the scale factor is given by a(t) ~ a(t s — t) . When t — > t 8 , therefore, 
a — > 00, namely, the Big Rip singularity appears. In this limit, the effective EoS 
parameter is given by 



1 16 



(26) 
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This behavior is identical with the case in which there is a phantom matter with 
its EoS being smaller than —1. As a consequence, we have reconstructed an explicit 
model with realizing a crossing of the phantom divide. 

From Eq. (fT9|) . we see that the effective EoS w e g becomes —1 when H = 0. By 
solving w c s — — 1 with respect to t by using Eq. (l20l) . namely, U(t) = 0, we find 
that the effective EoS parameter crosses the phantom divide at t = t c , given by 

1/(27+1) 



t c = t s [ -2 7 +./4 7 2 + 



7 



(27) 



It follows from Eq. (|2lj) that when t < t c 
time derivative of U(t) is given by 

,2 



7+1, 

U(t) > because 7 > 0. Moreover, the 



dU(t) 
dt 



2 7 ( 7 +l)(2 7 +ir 



15 



7+ ft + !)(£) 



27+1 



27 



27+1 



(28) 



Eq. (|28[) implies that the relation dU(t)/(dt) < is always satisfied because we 
consider only the period < t < t s as mentioned above. This means that U(t) 
decreases monotonously. Thus, the value of U (t) evolves from positive to negative. 
From Eq. ([20|) . we see that the value of w e g crosses —1. Once the universe enters the 
phantom phase, it stays in this phase, namely, the value of w e g remains less than — 1, 
and finally the Big Rip singularity appears because U(t) decreases monotonically. 
We note that there could be other types of the finite-time future singularities in 
modified gravity as shown in Refs. |33j |34] 

By using Eqs. ([13]). ([13]) . (fT5|) and <jTTJ) , P{t) is obtained as 

(4)' 15 



P{t) = 



It follows from Eqs. (TT2"|) and 



1 



\ 27+1 



that Q(t) is given by 



3=± 



(29) 



Q(t) 



-6H 



1 



(4) 



27+1 



E 

j=± 



-H 



t 



(30) 



If Eq. (|2^|) can be solved with respect to t as t — t(R), in principle we can find the 
form of F(R) by using this solution and Eqs. ((H), (|2"9"|) and (|3"0")) . However, for the 
general case it is difficult to solve Eq. (|22l) as t = £(i?). Therefore, as an solvable 
example, we illustrate the behavior of tlF(R) as a function of R = t^i? in Fig. 1 for 
7 = 1/2, p+ = -l/tf + , p_ = 0, /3+ = (l + 2 V / 19) /2 and t s = 2t . The quantities 
in Fig. 1 arc described in dimensionless quantities. The horizontal and vertical axes 
show R and t 2 s F, respectively. (Here, R = t 2 s R — 4R/Ro, where i?o is the current 
curvature. In deriving this relation, we have used t s = 2< , to s» H^ 1 , where H is 
the present Hubble parameter.) From Fig. 1, we see that the value of F(R) increases 
as that of R becomes larger. 



March 3, 2010 10:39 WSPC/INSTRUCTION FILE Bamba-6 



Crossing of Phantom Divide in F(R) Gravity 7 



800 



600 



% 400 



200 



200 400 600 800 1000 
R 



Fig. 1. Behavior of t 2 s F(R) as a function of R for 7 = 1/2, p + = -l/i'f + , p- = 0, /3+ 
(l + 2 v / l9) /2 and t s = 2t . 



To explore the analytic form of F(R) for the general case, we examine the 
behavior of F(R) in the limits t — > and t — > t s . When t — > 0, from Eq. (j23|) we get 



607(207- 1) 



7? 



(31) 



In this limit, it follows from Eqs. (g]), ^ and ^ that the form of 

F(R) is given by 



F(R)r 



^607(207 -l)^- 1 ^ 



to 



1 - 



[^V607(207- 1R 



1/2 



27+1 



57-1-/3 



2O7- 



pj [6O7 (20 7 - l)] ft/2 i?-^ /2 | . (32) 



Note that such action belongs to general class of actions with positive and neg- 
ative powers of curvature introduced in Ref. 1351 

On the other hand, when t — > t s , from Eq. (|25|) we obtain 
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In this limit, it follows from Eqs. Q, (gH), §U§ and §5§ that the form of 

F(R) is given by 

1 D Y.pA ts ~ 3v ^ R 



1 - 



1 _<wm R -i/2 



2 7 +l 



R 



ft' 



J 1 / 20 

x r-VT 



-1/2 



1 



f s - 3VT40R- 1/2 



(34) 



For large i?, namely, t 2 R^> 1, the expression of ^(i?) in (|34[) can be approximately 
written as 



F(R) « - 



1 



3V140(27+ 1) \*o 



(35) 



3. Summary 

We have studied a crossing of the phantom divide in F(R) gravity. We have re- 
constructed an explicit model of F(R) gravity in whic h a cro ssing of the phantom 
divide can occur by using the reconstruction methodP^EH As a result, we have 
shown that the Big Rip singularity may appear in the reconstructed model of F(R) 
gravity. We finally mention that by adding R 2 term (as it was first proposed in 
Ref. [24)) to the model or by adding non-sing ular theory^ R 2 (R n + Ci) / (R n + c 2 ), 
where n, c\ and ci are constants, B ig Rip singularities could be avoided because 
phantom behavior becomes transient pSCSI 
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